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ABSTRACT 

A high intensity source of a single neutrino flavour with known spectrum is most desirable 
for precision measurements, the consensus direction for the future. The beta beam is an 
especially suitable option for this. We discuss the prospects of a very long baseline beta 
beam experiment with a magnetized iron calorimeter detector. In particular, with the source 
at CERN and the detector at the proposed India-based Neutrino Observatory (INO) the 
baseline is near the 'magic' value where the effect of the CP phase is small. We observe that 
this experiment will be well suited to determine the sign of m\ — m| and will be capable of 
probing 9 13 down to about 1°. 
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I Introduction 



There is now compelling evidence in support of neutrino mass and mixing pQ from a number 
of atmospheric |2|, solar [Hj, reactor and long-baseline |Hj neutrino experiments. The 
neutrino mass eigenvalues and the Pontecorvo, Maki, Nakagawa, Sakata (PMNS) mixing 
matrix [S| provide a natural framework for formulating the scenario for three active neutrinos. 
At present, information is available on two neutrino mass-square differences and two mixing 
angles: From atmospheric neutrinos one gets the best-fit values with 3a error 1 |Am| 3 | ~ 
2.12lJ;g? x 1(T 3 eV 2 , 9 23 ^ 45.0 o ^ 3 p o while solar neutrinos tell us Am 2 12 ~ 7.9 x 1(T 5 eV 2 , 
9 12 — 33.21° [7j. At the moment, the sign of Am| 3 is not known. It determines whether 
the neutrino mass spectrum is direct or inverted hierarchical. The two large mixing angles 
and the relative oscillation frequencies could be useful for measurement of CP-violation in 
the neutrino sector, if the third mixing angle, #13, and the CP phase, 5, are not vanishingly 
small. The current bound on the former is sin 2 9i3 < 0.05 (3er) [HI Ej while the latter is 
unconstrained. 

A number of possible high-precision neutrino oscillation experiments are being designed to 
shed light on #13, 5, and the sign of Am| 3 : Among these are super-beams (very intense 
conventional neutrino beams) [TQ1 EJ 5 neutrino factories (neutrino beams from boosted- 
muon decays) ^3] , improved reactor experiments [Hj , and more recently (3 beams (neutrinos 
from boosted-ion decays) [To"! ITol IT7|. 

Here we focus on a long baseline (~ several thousand km) (3 beam experiment in conjunc- 
tion with a magnetized iron calorimeter detector with charge identification capability. The 
proposal for a detector of this type (ICAL) is being evaluated by the INO collaboration 
|18j . We consider the beta beam source to be located at CERN. To maintain collimation 
over such long baselines, the beta beam has to be boosted to high 7. The longer baseline 
captures a matter-induced contribution to the neutrino parameters, essential for probing the 
sign of Am^. The CERN-INO distance happens to be near the so-called 'magic' baseline 
[T9*] for which the results are relatively insensitive to the yet unconstrained CP phase. This 
permits such an experiment to make precise measurements of the mixing angle #13 avoiding 
the degeneracy issues [2D] which plague other baselines. 

II The (3 beam 

The beta beam, an idea put forward by Zucchelli [H], is connected with the production of a 
pure, intense, collimated beam of electron neutrinos or their antiparticles via the beta decay 
of accelerated radioactive ions circulating in a storage ring |21j . In particular, such a beam 
can be produced with the help of the existing facilities at CERN. An intense proton driver 
and a hippodrome-shaped decay ring are the essential requirements for this programme. 
It has been proposed to produce v e beams through the decay of highly accelerated 18 Ne ions 
and v e from 6 He [21]. Using the SPS accelerator at CERN, it will be possible to access 7 
~ 100 for completely ionized 18 Ne and 7 ~ 60 for 6 He. The ratio between the two boost 

1 Here Am?, ■ = m? — m? . 
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factors is fixed by the necessity of using the same ring for both ions. It is envisaged to 
have both beams simultaneously in the ring. Such an arrangement will result in a v e as well 
as a v e beam pointing towards a distant target. Higher values of 7, as required for longer 
baselines to INO from CERN, for example, can be achieved by upgrading the SPS with 
superconducting magnets or by making use of the LHC The reach of the LHC will be 7 = 
2488 ( 6 He) and 7 = 4158 ( 18 iVe) [221 • The beta beam is almost systematic free. 
u e are produced by the super-allowed (3~ transition \He —ALi + e~ + v e . The half-life of 
\He ++ is 0.807s and the Q value of the reaction is E Q =3.507 MeV. Neutrino beams can be 
produced by the super-allowed f3 + transition ioiVe -^^F + e + + v e , having the half-life 1.672s 
and the Q value, E =3.424 MeV. According to feasibility studies [23 Ell, the number of 
injected ions in case of anti-neutrinos can be 2.9 x 10 18 /year and for neutrinos 1.1 x 10 18 /year. 



Ill Neutrino fluxes 



Neglecting small Coulomb corrections, the differential width of /5-decay is described by: 
d 2 T* 1 In 2 
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where m e is the electron mass and E* is the neutrino energy 2 . Here E represents the electron 
end-point energy, tyz is the half life of the decaying ion in its rest frame and 
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where y e = m e /E . 

Since the spin of the parent nucleus is zero, it decays isotropically and there is no angular 
dependence in its rest frame. The Jacobian, J = [7(1 — /3cos0)] _1 , connects the rest frame 
quantities (cos9*,E*) to the lab frame ones (cos 9, E u ). 
The flux N is related to T by the radioactive decay law 



d 2 N 
dE v dt 



gjr 



dT 
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where g is the number of injected ions per unit time and r is the lifetime of that ion in its 
rest frame. 

We replace dVl by ^4, where dA is the small area of the detector and L is the distance between 
the source and the detector. So, using eqs. [T] and El the number of electron neutrinos, within 
the energy range E v to E v + dE u , hitting unit area of the detector located at a distance L 
aligned with the straight sections of the storage ring in time dt is given by: 
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where E* = jE u (l - [3 cos 6). 

From a technical point of view, it is not difficult to achieve designs aiming at higher 7 by 
direct extrapolation of existing facilities [22l 123] . The neutrino parameters we are interested 
to explore require a "high" 7 option (7 > 1500) which would be accessible, as noted earlier, 
in the LHC era at CERN. 

We discuss below the physics reach of a /3 beam using a magnetized iron calorimeter detector, 
of the type being considered by the India-based Neutrino Observatory (INO) collaboration. 
For the long baselines suitable for a rich physics harvest, the iron calorimeter detector (ICAL) 
being examined for INO provides a favourable target. The site for this Observatory has been 
narrowed down to one of two possible locations (a) Rammam in the Darjeeling Himalayas 
(latitude = 27°2'N, longitude = 88°16'E) or (b) Singara (PUSHEP) in the Nilgiris (latitude = 
11°5'N, longitude = 76°6'E). In the following, we show that ICAL will be an attractive choice 
for a very long baseline (3 beam experiment, with the source at CERN, Geneva (L = 6937 
km (Rammam), 7177 km (PUSHEP)). The unoscillated neutrino and anti-neutrino fluxes 
reaching the detector are depicted in Fig. [T] 

It is noteworthy that the large CERN to INO distance ensures a significant matter effect 
contribution enabling a determination of the mass hierarchy. At the same time, it matches 
the so-called 'magic' baseline ^H] where the results become insensitive to the unknown CP 
phase 5. This permits a clean measurement of # 13 . 




1.2e+07 



10 12 14 16 18 20 
Energy (GeV) 




~% 6~ 
Energy (GeV) 



Figure 1: Boosted spectrum of neutrinos and anti-neutrinos at the far detector assuming no oscillation. 
The flux is given in units of yr _1 m _2 MeV _1 . 
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IV Three flavour oscillations 



Here we briefly summarize the notations and conventions that will be followed. The neutrino 
flavour states \u a ) (a = e, fi,r) are linear superpositions of the mass eigenstates \vi) (i = 1, 
2, 3) with masses rrii, i.e., \u a ) = ^ U ai \ui). Here U is the 3x3 unitary matrix parametrized 
as (ignoring Major ana phases): 

where 



U = V 23 W 13 V 12 , (5) 




C12 S12 0\ / ci 3 s 13 e- l) 

V 12 = | -S12 Cia , W 13 = 1 I ,V 23 = | C 23 S23 | ■ (6) 

1/ \s 13 e iS ci 3 

where c 12 = 0056*12, S12 = sin6 l 1 2 etc., and 5 denotes the CP-violating (Dirac) phase. The 
probability that an initial Vf of energy E gets converted to a v g after traveling a distance L 
in vacuum is 

P{ Vf - i/,) = 5 /9 - 4 ^ Re(U* fi U 9i U fj U* gj ) sin 2 (1.27Am^) 

± 2^Im(C/;.C/ fli C/ /i C/;.)sin(2.54Amj|) (7) 

In the above, L is expressed in km, E in GeV and Am 2 in eV 2 . The - (+) refers to neutrinos 
(anti-neutrinos) . 

Neutrino interactions in matter modify the oscillation probability. Interactions of the u e 
occur through both charged and neutral weak currents making an additional contribution 
to its mass while the muon- and tau-neutrinos get contributions only through the neutral 
interaction. This alters both the mass splittings as well as the mixing angles. The general 
expression for the oscillation probability is messy. The appearance probability {y e —>■ v^) in 
matter, upto second order in the small parameters a = Am 2 2 /Am 2 3 and sin26 l 1 3, is |25,: 

2 „ 2 „ sin 2 [(l-i)A] 
P eM ~ sin 2 2013 sin 2 6 23 - 1 1 



± a sin 2^i3^ sin 5 sin( A) 



+ a sin 2#i3 £ cos 5 cos(A) 



I- A) 2 

sin(iA) sin[(l - A) A] 



A (l-A) 
sin(iA) sin[(l - A) A] 
~~A (l-A) 



+ a 2 cos 2 023 sin 2 20i 2 ^-^^; (8) 

where A = Am 2 n L/(4E), £ = cos6 13 sin20i 2 sin20 2 3, and A = ±{2\/2G F n e E) / Am 2 13 . G F 
and n e are the Fermi coupling constant and the electron density in matter, respectively. The 
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sign of the second term is positive (negative) for u e — > [y^ — > u e ). The sign of A is positive 
(negative) for neutrinos (anti-neutrinos) with normal hierarchy and it is opposite for inverted 
hierarchy. We have checked numerically that for low 9\3 (^4°) the results from the above 
approximate expression agree well with those from the exact three flavor oscillation formula. 
For higher values of #13 though agreement of a qualitative nature remains, the actual results 
differ by upto ~ 35%. 

One of the complications which needs to be addressed in the extraction of the neutrino 
properties is the issue of parameter degeneracies [20J; namely, that different sets of values 
of these parameters can result in the same predictions. It is imperative therefore to identify 
situations where this degeneracy problem can be circumvented or evaded. For example, in eq. 
ISl if one chooses sin(AA) = 0, the 5 dependence disappears and thus a clean measurement 
of the hierarchy and #13 is possible without any correlation with the CP phase S [F3\. The 
first non-trivial solution for this condition is \[2GFn e L = 2n. For an approximately isoscalar 
(one electron per two nucleons) medium of constant density p this equation gives an estimate 
of the size of this 'magic' baseline L magic : 

L magic [km] « 32726 - (9) 
p[gm/cm J 

In particular, for the CERN-INO path, the neutrino beam passes through the mantle of the 
earth where the density can be considered to be constant to a reasonable accuracy. The 
appropriately averaged density turns out to be p = 4.15 gm/cc for which L magic = 7886 km. 
The results presented in the course of our discussion are obtained by numerically solving the 
full three-flavour neutrino propagation equation based on the framework of Barger et ai, 
|26j . including the CP phase 5 and reflect the expectations for a near-'magic' baseline. 
Simulation for the ICAL design has shown excellent energy determination and charge iden- 
tification capability for muons with the energies relevant here. We focus therefore on the 
muon neutrino appearance mode, i.e., u e — > and v e — > transitions. Even though it 
is possible to increase the ion energy to achieve the threshold necessary for r production 
[y T appearance), it would require a very large storage ring and an enhanced storage time 
because of the lifetime dilatation. 



V Cross sections, Detector 

Following the standard approach, the neutrino-nucleus interaction cross section is obtained 
by including contributions from the exclusive channels of lower multiplicity (quasi-elastic 
scattering |2Z| and single-pion production 28 ), while all additional channels are incorporated 
as part of the deep-inelastic scattering [2H! cross section: 

0"CC = OQE + CItt + (?DIS- (10) 

At the low energy end, quasi-elastic events are dominant and the cross section grows rapidly 
for E v < 1 GeV, while at the higher energies (E v > a few GeV), mostly deep-inelastic scat- 
tering occurs and the cross section increases linearly with neutrino energy. At intermediate 
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energies, both types of events contribute. In addition, resonant channels dominated by the 
A(1232) resonance [2E] also take part in the process. We include all of the above. Because 
the neutrino energy extends up to about 20 GeV most events are deep-inelastic. There is 
about 10% contribution of quasi-elastic and single-pion production events each. 
The detector is assumed to be made of magnetized iron slabs with interleaved active detector 
elements as in the MINOS [H], and proposed ICAL detector at INO [H]. For ICAL, glass 
resistive plate chambers have been chosen as the active elements. In these proposals the 
detector mass is almost entirely (> 98%) due to its iron content. Here we follow the present 
ICAL design - a 32 Kt iron detector with an energy threshold around 800 MeV. The signature 
for the v e — > and v e — > transitions is the appearance of prompt muons whose tracks 
inside the detector will be reconstructed to give the direction and energy. Simulations have 
shown that the charge identification efficiency is around 95% at ICAL. So and will 
be readily distinguished. For this analysis the detector is taken to be of perfect efficiency 
and with no backgrounds 3 . Since we are very far from the source and the storage ring, the 
geometry of the storage ring will not play a vital role in the calculation of flux. 

VI Results 

VI. 1 Determination of the sign(Ara2 3 ) 

The CERN to INO distance is close to the 'magic' baseline (~ 7000 km) where matter effects 
are significant and the impact of the CP phase is negligible 4 . Over such long baselines, 
measurement of the neutrino mass hierarchy becomes possible, as matter effects become 
sizable. Within the three neutrino mixing framework, the results on solar neutrinos prefer 
the dominant mass eigenstates in v e to have the hierarchy m 2 > mi so that the mass- 
squared difference Am^ = m\ — m\ > 0. The sign of Am^ is a remaining missing piece 
of information to pin down the structure of the neutrino mass matrix. The beta beam can 
make good progress in this direction. 

Fig. |21 shows the number of events over a five-year period as a function of #13, taking the 
direct m| > 0) and inverted {m\—m\ < 0) hierarchies 5 . The noteworthy feature of this 
analysis is that for the direct hierarchy, the number of events obtained from a neutrino beam 
could be substantial while that for the anti-neutrino beam would be strongly suppressed, 
while the opposite will be true for the inverted hierarchy. Such an asymmetry would be easy 
to detect using the charge identification abilities of ICAL. 

The mass hierarchy can be probed at the 4.4 (4.8)<r level with a neutrino (anti-neutrino) beam 
for values of #13 as low as ~ 1°. As seen from fig. the sensitivity increases dramatically 
with #13. This sensitivity will also depend on the precise value of Am|. For example, for 

3 Atmospheric neutrino and other backgrounds will be eliminated by the directionality cut imposed in 
event selection. 

4 For the results which we present, we have used the CERN to Rammam (L = 6937 km) baseline. We 
have checked that if the baseline for the alternate PUSHEP site (L — 7177 km) is used, the results vary by 
less than 5%. 

5 Unless specified otherwise, where necessary, we use the best-fit values of the mixing parameters. 
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Am 2 ^ within the present la interval [1.85 - 2.48] xl0~ 3 eV 2 , this significance varies within 
3.5 - 5.3cr (4.6 - 5.1a) for neutrinos (anti-neutrinos). In the above, the CP phase 5 is chosen 
to be 90°. As checked below, the CERN-INO baseline, close to the 'magic' value, ensures 
essentially no dependence of the final results on 5. 

In this calculation, we have considered an uncertainty of 2% [TH] in the knowledge of the 
number of ions in the storage ring. Following the standard practice, we have assumed a 
10% fluctuation in the cross section, a. The statistical error has been added to the above in 
quadrature. We have neglected nuclear effects. 
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Figure 2: The number of events as a function of #13 for neutrinos (anti-neutrinos) is shown in the left 
(right) panel. The solid (broken) curves correspond to Am| 3 < (Am| > 0). 



VI. 2 Precision measurement of # 13 

Aside from the neutrino mass hierarchy, the other major unknown in the neutrino sector 
is the mixing angle 6 13 . Here also the results for the long baseline beta beam set-up are 
encouraging. 

In Fig. El we plot the number of events in 5 years as a function of #13 for two extreme values 
of S. Results are shown for neutrinos and anti-neutrinos. The dependence on 5 is seen to 
be very mild - a reflection of the near 'magic' baseline. The growth of the number of events 
with increasing # 13 is consistent with eq. (JHJ). For these plots we have chosen Am^ > 0. In 
this already noted in Fig. EJ the number of events for the v e beam is quite small 

while for v e it is substantial. Therefore, for this mass hierarchy, the neutrino run must be 
used to extract 6*13 and values as small as 1° can be probed. For the opposite hierarchy, the 
anti-neutrino beam will give the larger number of events which can be used to determine 

013- 
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The estimated 3a errors on #i 3 measured to be 1°(5°) are tgjo (ti^o) with 5 = 0° for neutri- 
nos. The results are somewhat worse for anti-neutrinos for the direct hierarchy, Am^ > 0, 
considered here. For the inverted hierarchy, anti-neutrino beams provide the better mea- 
surement. 

In the extraction of #13, a major role is played by the value of Am^. For illustration, with 
a neutrino beam and for 5 = 90°, the la error of Am^ translates to uncertainties of ~ ±1° 
at #13 = 5° and less than ±| at #13 = 1°. This is for the normal hierarchy. In principle, 
the long baseline beta beam experiment can narrow down the permitted range of Am^. 
However, it is very likely that this improvement will be achieved in the meanwhile by other 
experiments. 

The effect of 5 is negligibly small for neutrinos as seen in Fig. |3] (left). To estimate the 
effect of 5 for the case of anti-neutrinos we vary it over its whole range and find that the 
uncertainty range is less than 1° for all #13. 




Figure 3: Variation of the number of events with #13 for v (left) and v (right) for a 5-year run. Here, 
Am,23 ' s cnosen positive. 

VII Conclusions 

We have discussed the prospects of obtaining information on the mixing angle (9 13 and the 
sign of Am^ using a magnetized iron calorimeter detector, such as the proposed ICAL 
detector at INO, and a high 7 beta beam source. It appears that such a combination of a 
high intensity z/ e , v e source and a magnetized iron detector is well-suited for this purpose. 
We have focused on the CERN to INO baseline, which is close to the 'magic' value, and 
found that it should be possible to determine the sign of Am^ and probe #13 down to 1° in 
a five-year run. The effect of the CP phase 5 is quite mild. 
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